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Localized structures in broad area VCSELs 
experiments and delay-induced motion 


Mustapha Tlidi, Etienne Averlant, Andrei Vladimirov, Alexander Pimenov, 
Svetlana Gurevich and Krassimir Panayotov 


Abstract We investigate the space-time dynamics of a Vertical-Cavity Surface- 
Emitting Laser (VCSEL) subject to optical injection and to delay feedback control. 
Apart from their technological advantages, broad area VCSELs allow creating lo¬ 
calized light structures (LSs). Such LSs, often called Cavity Solitons, have been 
proposed to be used in information processing, device characterization, and others. 
After a brief description of the experimental setup, we present experimental evi¬ 
dence of stationary LSs. We then theoretically describe this system using a mean 
field model. We perform a real order parameter description close to the nascent 
bistability and close to large wavelength pattern forming regime. We theoretically 
characterize the LS snaking bifurcation diagram in this framework. The main body 
of this chapter is devoted to theoretical investigations on the time-delayed feedback 
control of LSs in VCSELs. The feedback induces a spontaneous motion of the LSs, 
which we characterize by computing the velocity and the threshold associated with 
such motion. In the nascent bistability regime, the motion threshold and the velocity 
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of moving LSs depend only on the feedback parameters. However, when consider¬ 
ing the previously introduced mean-field model, theoretical predictions indicate that 
both motion threshold and velocity are strongly affected by the phase of the delay 
and by the carrier relaxation rate. 


1 Introduction 

Spontaneous symmetry breaking and self-organization phenomena have been ob¬ 
served in various fields of nonlinear science such as nonlinear optics, fibre optics, 
fluid mechanics, granular matter or plant ecology. The link between the well known 
Turing instability and transverse patterns formation in nonlinear optics was estab¬ 
lished for the first time by Lugiato and Lefever ||T|. In their seminal paper, they 
considered an optical resonator filled with a passive nonlinear medium and driven 
by a coherent radiation beam. Since then, many driven systems have proven to allow 
periodic patterns in the transverse section of their output beam. Besides a periodic 
distribution of light, LSs may form in the plane perpendicular to the propagation 
axis. They are often called localized spots and localized patterns, or cavity soli- 
tons which appear either isolated, randomly distributed or self-organized in clus¬ 
ters forming a well-defined spatial pattern EH. When LSs are sufficiently sep¬ 
arated from each other, localized peaks are independent and randomly distributed 
in space. However, when the distance between peaks decreases they start to inter¬ 
act via their oscillating, exponentially decaying tails 01216117]. LSs have been 
reported in nonlinear resonators such as lasers with saturable absorbers Emm, 
in passive nonlinear resonators El 12 HU, optical parametric oscillators 112 in 
left-handed materials iH [12 ns El, in exciton-polariton patterns in semiconduc¬ 
tor microcavities lITSl E] and in the framework of the Ginzburg-Landau equation 
EqI EH 122 E2 iMl- Phase solitons have been demonstrated far from any pattern 
forming instability ll25l l26l l27l l28l l29l [30l . 

Localized structures are homoclinic solutions (solitary or stationary pulses) of 
partial differential equations. The conditions under which LSs and periodic patterns 
appear are closely related. Typically, when the Turing instability becomes subcrit- 
ical, there exists a pinning domain where LSs are stable. This is a universal phe¬ 
nomenon and a well-documented issue in various fields of nonlinear science. The 
experimental observation of LSs in driven nonlinear optical cavities has further mo¬ 
tivated the interest in this field of research. In particular, LSs could be used as bits 
for information storage and processing. Several overviews have been published on 
this active area of research iCT 1^ 1^ [3213113213^ l38l [3^ 1401 l4Tl 1^ l43l l44l 
|45l|46ll42|48|. 

Many theoretical and experimental studies on LS formation in VCSELs have 
been realized Il42l50l . They have been experimentally observed in broad area VC¬ 
SELs both below li49] [STIl and above lf52l the lasing threshold when injecting a 
holding beam with appropriate frequency and power. A spatially LS has also been 
found in a medium size VCSEL, but only by using its particular polarization proper- 
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ties 1531. Cavity soliton lasers (CSLs) in a VCSEL system without a holding beam 
have been demonstrated both experimentally l54l and theoretically 1551 in VCSELs 
subject to frequency selective optical feedback and in face to face coupled VCSELs 
l5^l57l . In these systems, the VCSELs are placed in self imaging optical systems 
with either an external grating or another VCSEL biased below lasing threshold, so 
that the system becomes bistable. Lasing spots spontaneously appear in these sys¬ 
tems and can be switched on and off by another laser beam. As a matter of fact a 
broad area VCSEL with saturable absorber has been the first system in which LSs 
have been predicted and studied theoretically l58ll59ll60l . LSs in a monolithic opti¬ 
cally pumped VCSEL with a saturable absorber have been demonstrated in ED and 
their switching dynamics studied in 1621 . Several applications of LSs in VCSELs 
have been demonstrated; optical memory E3l . optical delay line 1^ and optical 
microscopy l50l . 

In this chapter, we investigate the formation of LSs in Vertical-Cavity Surface- 
Emitting Lasers (VCSELs) subject to both optical injection and delay feedback con¬ 
trol. These lasers are characterized by a large Lresnel number and a short cavity. 
VCSELs are the best candidate for localized structures formation in their transverse 
section. The first VCSELs were fabricated in 1979 1651 and later on they reached 
performances comparable to those of edge-emitting lasers l66llC7l . Nowadays VC¬ 
SELs are replacing edge-emitting lasers in short and medium distance optical com¬ 
munication links thanks to their inherent advantages; much smaller dimensions, cir¬ 
cular beam shape that facilitates coupling to optical fibres, two-dimensional array 
integration and on wafer testing that brings down the production cost EtI . As VC¬ 
SELs emit light perpendicular to the surface and the active quantum wells, their 
cavity length is of the order of 1 /rm - the wavelength of the generated light. Thanks 
to the maturity of the semiconductor technology VCSELs can be made homoge¬ 
neous over a size of hundreds of /rms while the characteristic LS size is about 10 
jum. The timescales of the semiconductor laser dynamics and LS formation are in 
the ns scale, which allows for fast and accurate gathering of data. Linally, VCSEL 
physics and dynamics are quite well understood llC7ll^l69]| . 

We investigate experimentally and theoretically the formation of stationary LSs 
in VCSELs and we describe theoretically the effect of a time delayed feedback 
on the stability of LSs. Lor this purpose, we adopt the Rosanov-Lang-Kobayashi 
approach for modelling of delay feedback lEQlllD- 

This chapter is organized as follows. After an introduction, we provide a de¬ 
scription of the VCSELs in Sec (2), the experimental setup and the observation of 
localized structures in a medium size VCSEL are described in Sec. (3). A theoret¬ 
ical description based on mean held model and the derivation of the generalized 
Swift-Hohenberg equation are given in Sec. (4). Stationary LSs and their bifurca¬ 
tion diagram are presented in Sec. (5). LSs brought into motion under the effect of 
delay feedback are discussed in Sec. (6). We conclude and draw some perspectives 
of the present work in Sec. (7). 
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2 Vertical-Cavity Surface-Emitting Lasers 

The structure of a VCSEL is by far more complex than the one of an edge-emitting 
semiconductor laser. But this complexity, and the inherent fabrication costs, did not 
keep VCSELs from becoming the second most produced laser type GZllTlllTl. 
The biggest advantages of this structure are the circular emission pattern and low 
divergence, which allow to easily couple the light to an optical fibre for optical 
interconnects. Moreover, the high reflectivity of the Bragg mirrors provide a low 
lasing threshold, which makes VCSELs the choice lasers for low power applications 
(such as optical mice for example). For high power applications, VCSELs can easily 
be coupled into arrays ll75l . which provides a high power coherent low divergence 
beam. 

The emission surface of a VCSEL can be made from a few /rm to some 200 
jLim - the first ones are best suited for low power applications, whereas the second 
category is typically used for spectroscopy. Apart of their use in spectroscopy, broad 
area VCSELs are best suited for LSs studies. They have several characteristics that 
makes them a choice material: 

• Their use in information technologies provides a mastering of the fabrication 
process, leading to a la carte properties; 

• The response time is much faster than the one in liquid crystals or photorefrac- 
tive media, which makes them more likely to be used in potential high speed 
applications; 

• Experiments can be carried out on a single one hundred micrometer chip, in 
comparison with some meter needed for creating LSs in gas cells; 

The broad-area bottom emitting VCSEL structure we use in our experiments is de¬ 
scribed in Fig. [T] The top and bottom distributed Bragg reflectors consist of 30 and 
20.5 Alo.ssGao.nAs-GaAs layer pairs, respectively. The active region consists of 
three Ino, 2 Gao gAs quantum wells embedded in GaAs barriers and AlGaAs cladding 
layers ll75l . The bottom-emitting configuration allows the stand-alone VCSEL to 
have a better (more homogeneous) current distribution in the transverse plane and, 
therefore, is more suitable for producing LSs. Moreover a heat sink can be directly 
mounted on top of the p-contact and the p-doped DBR, which is the main source of 
heat. The temperature distribution is hence more homogeneous as well. 


3 Experimental observations 

3.1 Description of the experimental setup 

The experimental setup used for the generation of two-dimensional LSs is shown in 
Fig .|^ The injection part consists of a tunable External Cavity Diode Laser (Master) 
in a Littrow configuration, optically isolated from the rest of the experiment. The 
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Fig. 1 View of a bottom- 
emitting VCSEL. The blue 
and green layers constitute 
the pairs in the Bragg mirrors. 
The red element is the cavity 
containing the quantum wells. 
The grey part is the substrate 
on which the structure has 
been grown and through 
which the light is emitted. 


Fig. 2 Experimental setup 
schematic. The full line is the 
path of the light from the mas¬ 
ter laser, whereas the dashed 
line is the path followed by 
the light from the VCSEL. 

(i): injection preparation and 
monitoring; Master: master 
laser, OI: optical isolator, A/2 
: half wave plate, M: mirror, 
VODF: variable optical den¬ 
sity filter; (ii) : VCSEL; (iii): 
analysis branch; PD: photo¬ 
diode, OSA: optical spectrum 
analyser. Reprinted from (76) 



Top contact 
Inner cavity and 
ctive layers 


Bragg mirrors 




linear polarization of the injected light is tuned to match the one of the VCSEL. In 
order to produce LSs, both the optical injection power and the detuning between the 
injection frequency and the VCSEL resonance frequency are adjusted. 

The VCSEL we use is a 80/im diameter, bottom emitting InGaAs multiple quan¬ 
tum well VCSEL, as described in Sec.(|^ which has a threshold cutTent of 42.5 mA 
at25.0°C. 

We analyse the optical spectrum, near field and output power of the VCSEL. 

A photography of the actual setup is shown in Eig[^ 


3.2 Experimental observations 

Tuning the width of the injection beam to 100 pm, the detuning to — 174GHz, the 
VCSEL driving cutTent to 45.013mA and its substrate temperature to 25.01°C, we 
obtain the bistability curve depicted in Eig.|^ We continuously varied the optical in¬ 
jection power by increasing it, before decreasing it. The hysteresis phenomenon as¬ 
sociated with this experiment is evidenced in Eig|^). The insets represent near held 
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Fig. 3 Photography of the experimental setup. Left down side of the picture: injection prepara¬ 
tion with external cavity laser diode, optical isolator, half-wave plate and variable optical density 
filter. Top right: injection monitoring and VCSEL. Down right: analysis branch; optical spectrum 
analyser, CCD camera and photodiode. 


profiles on the higher and lower branch of the hysteresis curve. One dimensional 
prohles along the horizontal lines drawn on the insets are shown in Fig|^)(up) 
and c)(down). The lower branch corresponds to the ’’homogeneous” steady state. In 
Fig|^ the difference between these two states is evidenced. On the higher branch, 
there is not only a peak, but a clearly noticeable oscillating tail around this peak, 
which is characteristic of a LS. 



Fig. 4 Bistability between one and two-peaked LSs inside the near field of the VCSEL as a func¬ 
tion of the optical injection power, (a): power emitted by the VCSEL as a function of the optical 
injection power for 6 = — 174GHz and a beam waist of 100pm. The insets (i) and (ii) respectively 
represent near field profiles on the higher and lower branch of the hysteresis curve, (b) and (c): one 
dimensional profiles along the horizontal line drawn on the aforementioned insets. Redrawn from 
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Fig. 5 Cross sections along 
the solid lines indicated in 
Figj^a), (i) and (ii). The 
dashed line is the state 
(ii)(lower branch of the hys¬ 
teresis), whereas the full line 
is the system with a LS (up¬ 
per branch of the hysteresis). 
Redrawn from CH 



4 Mean field model 


In this section, we describe the dynamics of a VCSEL submitted to optical injec¬ 
tion and delayed feedback. To do so, we assume an external cavity (i) in which the 
diffraction is compensated (ii) much longer than the characteristic diffraction length 
of the field. We further apply (i) a paraxial approximation (ii) a slowly varying en¬ 
velope approximation. The mean field model describing the space-time evolution 
of the electric field envelope E and the carrier density in a VCSEL subjected to 
optical injection and time-delayed feedback is then given by the following set of 
dimensionless partial differential equations 


dE 

Ii 


{lJ.+i9)E + 2C{l-ia){N-l)E 


dN 

dt 


+ Ei- 
= -7 


^e'‘^E{t-x)+N^E, 

N -1 + {N -\)\E\^ - dy^N 


( 1 ) 

( 2 ) 


The parameter a describes the linewidth enhancement factor, /i and 0 are the cavity 
decay rate and the cavity detuning parameter, respectively. The parameter £; is the 
amplitude of the injected field which we assume to be positive in order to fix the 
origin of the phase. C is the bistability parameter, y is the carrier decay rate, I is 
the injection current, and d is the carrier diffusion coefficient. The light diffraction 
and the carrier diffusion are described by the terms N'^E and dV^N, respectively, 
where is the Laplace operator acting in the transverse plane (x, y). Below we 
consider the case when the laser is subjected to coherent delayed feedback from 
an external mirror. To minimize the effect of diffraction on the feedback field we 
assume that the external cavity is self-imaging CD. The feedback is characterized 
by the delay time T = 2Lext / c, the feedback rate rj >0, and phase (p, where Lgxt is the 
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external cavity length, and c is the speed of light. The link between dimensionless 
and physical parameters is provided in ll78l . Using the expression for the feedback 
rate rj = r^/^(l —/?)//?given in ESI, where r (R) is the power reflectivity of 
the feedback (VCSEL top) mirror and T,„ is the VCSEL cavity round trip time, we 
see that the necessary condition for the appearance of the soliton drift instability 

Rt 

Tj T > 1 ifTTl can be rewritten in the form r > 

T = 20t,>, the latter inequality becomes r > 1.5 


T. In particular, for R = Q3 and 


10 -^ 

To reduce the number of parameters, we introduce the followi ng ch ange of vari- 

^ of a VCSEL 


ables; n = [2C{N — 1) — l]/2 and e — E*/\/2. The model Eqs. (1 
driven by an injected field Y = Eil{2\/2) take the following form: 


dte = W'e + {I + ia)ne+ Y + rj'e e{t — x) — (3) 

(9;n = 7[P —n—(1+2n)|ep+DV^n]. (4) 


The pump parameter P is P = C{I — 1) — 1/2, y = "/12, D — 2d, rj' = ^/2, and 
9' = (9 + a)/2. The new time and space scales are (f,T) = 2{t',x') and = 2V'^. 
Let us assume for simplicity that the detuning is 0' = 0 and the feedback phases are 


\jf = 0 or xjf ^ 7t. 

The homogeneous steady states of Eqs. (3 4i are T = —es{l+i(x){P— \es'\^)/[l + 
2|ejp) and = {P — \es'\^)/{ \ +2\es\^). We explore the vicinity of the nascent 
optical bistability regime where there exists a second order critical point marking 
the onset of a hysteresis loop. The critical point associated with bistabilty is obtained 
when the output intensity as a function of the injection parameter Y has an infinite 
slope, i.e., dY/d\es\ = d'^Y/d\es\^ = 0. The coordinates of the critical point are 
ec = {1- /a)y/3/2(l + a2), = -3/2, P^ = -9/2, = 8a/[3(l + a^)] and 

T; = (3/2)(3/2)( 1 + Ct^) Our objective is to determine a slow time and slow space 
amplitude equation which is valid under the following approximations: (i) close to 
the onset of bistabilty (ii) close to large wavelength symmetry breaking instability. 
Starting from Eqs. Q |^, the deviation u of the electric field from its value at the 
onset of bistabilty is shown to obey | 


dfU = y — u{p + u^) + r\u(t — x) 

+ {d-^)y^u-aV'‘u-2{Vuf, 


(5) 


where a = (1 — a^)/(4a^). The parameter y denotes the deviation of the injected 
field amplitude from Yc- The real variable u, the parameters p and d are the devi¬ 
ations of the electric field, the pump parameter and the carrier diffusion coefficient 
from their values at the onset of the critical point, respectively. In the absence of 
delay; i.e., rj =0, Eq.(|^ is the generalized Swift-Hohenberg equation that has been 
derived for many far from equilibrium systems 118111821 [^ . Even in absence of the 
delay feedback term, the terms uV^u and (Vu)^ render Eq.lj^ nonvariational. 











Localized stmctures in broad area VCSELs: experiments and delay-indueed motion 


9 


5 Stationary localized structures 

LSs are nonlinear bright or dark peaks in spatially extended systems. Such structures 
have been observed in the transverse section of coherently driven optical cavities, 
and are often called cavity solitons. Currently they attract growing interest in optics 
due to potential applications for all-optical control of light, optical storage, and in¬ 
formation processing ma. When they are sufficiently separated from each other, 
localized peaks are independent and randomly distributed in space. However, when 
the distance between peaks decreases they start to interact via their oscillating, ex¬ 
ponentially decaying tails. This interaction then leads to the formation of clusters 
01211610. Mathematically speaking, LSs are homoclinic solutions (solitary or sta¬ 
tionary pulses) of partial differential equations. The conditions under which LSs 
and periodic patterns appear are closely related. Typically, when the Turing insta¬ 
bility becomes sub-critical, there exists a pinning domain where localized structures 
are stable. This is a universal phenomenon and a well documented issue in various 
fields of nonlinear science, such as chemistry, plant ecology, or optics (see some 
recent overviews on this multidisciplinary issue ( 141114211481 ). 

In this section we describe some basic properties of stationary LS and their bi¬ 
furcation diagrams in a one dimensional setting. In the absence of delay feedback, 
Eq. (|^ admits a variety of LSs. The generalized Swift-Hohenberg equation^ is nu¬ 
merically integrated using a classical spatial finite-difference method with forward 
temporal Euler integration. The boundary conditions are periodic in both transverse 
directions and the initial condition consists of a large amplitude peaks added to the 
unstable homogeneous steady state. We fix all the parameters except the amplitude 
of the injected field y. Examples of localized clusters having an odd or an even num¬ 
ber of peaks are shown in Eigs. lb and Ic. They are obtained for the same parameter 
values and differ only by the initial condition. 

We examine the case of one-dimensional monostable system in the subcritical 
regime where the homogeneous steady state coexists with a spatially periodic struc¬ 
ture. In addition, the system exhibits a high degree of multistability in a finite range 
of y values often called the pinning region ll84l . More precisely, Eq. (0 then ad¬ 
mits an infinite set of odd and even cavity solitons as shown in Eig. i.e., a 
set of stationary solutions that exhibit n — 2p — 1 or n = 2p peaks, where p is a 
positive integer. The limit p —> oo, corresponds to the infinitely extended periodic 
pattern distribution. In the pinning region, the width of LS is close to the half of 
the periodic wavelength structure. Since the amplitudes of localized patterns having 
different number of peaks are close to one another, in order to visualize the clusters 
properties, it is convenient to plot the “L 2 -norm” 

,yV = j dxdy\u — Us'^ (6) 

as a function of the injected field y. This yields the two snaking curves with odd 
or even number of peaks as shown in bifurcation diagram of Eig. ®). As in¬ 
creases, at each turning point where the slope becomes infinite, a pair of additional 
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peaks appears in the cluster. One sees that this behavior, referred to as homoclinic 
snaking phenomenon Il85ll8^[87ll88l[8^l90ll^ . and recently observed experimen¬ 
tally 192119^ . corresponds to back and forth oscillations inside the pinning region. 
In two spatial dimensions, the variety of stable localized patterns is much larger than 
in one dimensional system. 



Fig. 6 One dimensional LSs. (A) Multiple peaks localized structures with odd or even number of 
peaks obtained fory = —0.35. (B) Snaking bifurcation diagram showing two inter-weaved snaking 
curves: the branches (a), (b), (c), (d) and (e) correspond to states having 1, 2, 3, 4 and 5 peaks 
in fig. |§k. Redrawn from (80). The full and the broken lines correspond, respectively, to stable 
and unstable localized branches of LSs. The parameters are p = —0.7, d = —1.2, a = 0.75, and 
T = 7] = 0. 


6 Moving localized structures 

In this section we investigate the effect of a time-delayed feedback control on the 
stability of LSs in VCSELs. This delayed feedback loop compares the electric field 
at the current moment of time and its values at some time in the past. Recent 
studies that combined analytical and numerical analysis of the two-dimensional 
Swift-Hohenberg equation suggested that steady LSs can become mobile when 

T7T< 1 177]|g|80l. 

V=|^-(1-T7T), (7) 


I j^^uldxdy 
V u^dxdy 


with 
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Due to the rotational symmetry of the generalized Swift-Hohenberg equation Q, 
there is no preferred direction for the motion of a circularly symmetric localized 
structure. The instability leading to the spontaneous motion of a LS solution is a cir¬ 
cle pitchfork type of bifurcation. Therefore, the x axis can be chosen for the estima¬ 
tion of the velocity v. In this case we obtain M] = duo(r)/dx and 112 = d^uo{r) jdx^ 
where Mo(r) is the stationary localized structure solution. We recover the expression 
for the soliton velocity (0 that was obtained hrst in the case of the Swift-Hohenberg 
equation IfTTl . The speed formula 0 is valid for any localized pattern, regardless of 
the number of peaks it contains. The factor Q describes the spatial form of the local¬ 
ized pattern. This factor can be be calculated only numerically. In particular, for the 
parameter values y = —0.35, p = —0.7, d = —1.2, a = 0.75, we obtain Q = 1.44. 
The velocity Eq. Q divided by the factor Q is plotted as a function the delay time 
for a hxed value of the feedback strength as shown in Fig. The curve of the 
velocity has a maximum at T = 2/tj, which corresponds to the maximal velocity 
Vmax = Q'nl^- Note that the threshold rjT = I and the expression for the formula 
( 7 ]) have been found later on for the chemical reaction-diffusion type of equations 

|S|96l|92l. 

Numerical simulations of the Eq. 0 show indeed that single and many peaked 
LSs exhibit a spontaneous motion as predicted by the above theoretical analysis, as 
shown in Fig.[^ 

The velocity and the threshold rfi = 1 associated with the LS motion depend 
only on the delay parameters J] and T. This statement is valid only in nascent optical 
bistability regime,where the dynamics of the system are described by the real order 
parameter equation Q, and for a hxed feedback phase of 1 / = 0 or y/ = ;r. In what 
follows we will examine the combined role of the phase of the delay feedback \j/ 
and the carrier relaxation rate 7 in the framework of the full mean held model. 

Since our system is isotropic, the motion of localized structures occurs in an 
arbitrary direction. At the pitchfork bifurcation the stationary LS loses stability and 
a branch of moving LSs with the velocity v = |v| bifurcates from the stationary 
LS branch of solutions. The bifurcation point can be obtained from the hrst order 
expansion of the uniformly moving LS in power series of the small velocity v. Close 
to the bifurcation point, the uniformly moving LS can be expanded in power series 
in the small velocity v and through the solvability condition, we obtain the drift 
instability threshold ll^ |9^ 


with 


1 + 7 ’(^/c) 

y/T+W^ ^ cos [(p + arctan {a /c)] 

a = (rJ',V^2)-(rl,V^i). 

b = (V/I,r3), 

c = (W, ri) + (rl, ¥ 2 }, 


(8) 


(9) 


and 
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T 


Fig. 7 Velocity by unit of the factor Q of moving localized structures as a function of the time 
delay T for different values of the delayed feedback strength rj. 


W={Wl,V2,W3f = dxiXo,Yo,Nof (10) 

a translational neutral mode of the linear operator L, L^f = 0, while \j/^ = , v4 > v4 

is the corresponding solution of the homogeneous adjoint problem Xj/^ = 0. The 
real A'o(;ii:,y) and the imaginary lo(x,y) parts of the electric held Eq and the carrier 
density No{x,y) are the stationary axially symmetric LS prohles. They correspond 
to the time-independent solutions of Eqs. 0 and Q with T = 0. The coefficients 
a and b are calculated numerically using the relaxation method in two transverse 
dimensions. 
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From the expression of the threshold associated with the drift instability Eq. 
we see that the product rjT is not unity as in the case of the generalized Swift- 
Hohenberg equation, but depends strongly on the feedback phase (p and carrier 
relaxation rate 7 . We plot in Fig. the threshold t] associated with the drift in¬ 
stability as a function of the phase of the delay feedback. The numbers on top of 
the different curves are the different values of the carrier decay rate 7 . The car¬ 
rier relaxation rate strongly affects the threshold associated with the drift instabil¬ 
ity as shown in Fig. When we increase the carrier relaxation rate, the threshold 
associated with the moving LS gets higher. In addition, we see from Fig. 1 ^ that, 
whatever the value of the carrier relaxation rate, the drift instability occurs only 
within the subinterval {%,{„ — Kjl, + 7t/2) of the interval {(pmin — (Pmin + '^), 
where ()?„,■„ = — arctana is the delay feedback phase, corresponding to the lowest 
critical feedback rate Tjg”" = (1 -f 7 ^*h)/(T-\/l + 0 ^). Note that when tjq increases 
very rapidly when approaching the boundaries of these subintervals. In the limit of 
fast carrier response, 7 3^ 1 , and for zero feedback phase, ^ = 0 , we recover the ex¬ 
pression the threshold formula tjt = 1 that has been obtained in both variational 
Swift-Hohenberg equation Cll, and in a modihed nonvariational one II 8 OII . Note that 
at 7 —^ 00 , a ^ 0, and (p = — arctana the critical feedback rate appears to be smaller 

1 /O 

than that obtained for the real Swift-Hohenberg equation, tjt = (1 -f 


Fig. 8 Space time map of 
a single-peaked localized 
structure solution to eq. (|^. 
Parameters are p = —0.9, 
rf = -1.5,y = -0.5, 77 = 0.15 
and T = 15. Redrawn from 
Ref. Cool 



Space 


To estimate the velocity of moving localized structure, we expand the slowly 
moving localized solution in the limit of a small velocity v. The detailed calculations 
can be found in ll99l . The velocity v of LSs then obeys 


V = with Q = {l/x)^Jq/{r7]). 


( 11 ) 


The factor Q is important since it determines how fast the LS speed increases with 
the square root of the distance from the critical feedback rate. The Q factor depends 
on the delay feedback phase as shown Fig. 10 The 5rj denotes the deviation of the 
feedback strength from the bifurcation point associated with the moving LS. The 
coefficients r and qamq = asin^ + ccos (p, and r = /sin^-t-gcos with 
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9 


Fig. 9 Threshold associated with the drift bifurcation as a function of the feedback phase (p is 
plotted for different values of the carrier relaxation rate y. Parameters are ^ = 1.0, 0 = —2.0, 
C = 0.45, a = 5.0, t = 100, d = 0.052, £,■ = 0.8, 1 = 2. The values of the parameter y are shown 
in the figure. Redrawn from (99). 


f — {y^i^^xxx^o} (V ^2 ’ — ( 1/^3 7 ^7 (1/^2 7 ■ 

The coefficients a and c are defined in Eq. 



9 


Fig. 10 Q factor as a function of the feedback phase (p , for different values of the carrier relaxation 
rate 7 . The Q factor describes the growth rate of the LS velocity with the square root of the deviation 
from the critical feedback rate. Parameters are ^ = 1.0, 6 = —2.0, C = 0.45, a = 5.0, T = 100, 
d = 0.052, Ei = 0.8, 1 = 2. The values of the parameter 7 are shown in the figure. Redrawn from 

(99] 


Numerical simulations of the full model Eqs. and the generalized Swift- 
Hohenberg equation Eq. Q agree with the above theoretical predictions. Indeed, 
when we choose parameter values such as the system operates above the threshold 
associated with the motion of LS, a single-peaked LS exhibits a regular motion in 
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the transverse plane of the cavity as shown in Fig. [^(1 -dimensional setting, (Eq|^) 
and [^2-dimensional setting, Eqs Q and Q). 


time = 3600, 4000, 4400 




Fig. 11 Field intensity |£|^(top) and carrier density N (bottom) of a single peaked two-dimensional 
moving cavity soliton at different times. Parameters are C = 0.45, 6 = —2, a = 5, 7 = 0.05, 
d = 0.052, p = 2, Ei — 0.8, r = 200, t] = 0.07, (p = 3.5. Results obtained using Eqs.^andj^ 
Redrawn from EH). 


Note, however, that the motion of LS under the effects of delay feedback is not 
always regular. It has been shown recently that LSs can exhibit a temporal chaos: 
numerical simulations of a broad area VCSEL with saturable absorber subjected to 
time-delayed optical feedback have shown evidence of complex temporal dynamics 
of LSs 01011 . This spatio-temporal chaos is localized in space. More recently, in 
the absence of decay feedback, it has been shown that the VCSEL with saturable 
absorber without optical injection injection may exhibit LSs that drift and oscillate 
simultaneously, and a chaotic behavior 0102111030 . 


7 Conclusions and perspectives 

In this chapter we have investigated the formation of localized structures in a 
Vertical-Cavity Surface-Emitting Laser subject to optical injection. This device con¬ 
sists of a medium size bottom-emitting InGaAs multiple quantum well VCSEL 
operating in an injection locked regime. In this regime, we have described exper¬ 
imentally the formation of stationary localized structures of light in the transverse 
section of this device. The experimental part has been performed in an injection 
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locked regime and in the absence of delay feedback control. We have characterized 
LSs by drawing their bifurcation diagram and performed a numerical simulation for 
the full model and the generalized Swift-Hohenberg equation. 

Then we have described the space-time dynamics of a VCSEL by adding the 
delay feedback control in the modelling by adopting a mean-field approach. The 
time-delayed feedback is modelled following a Rosanov-Lang-Kobayashi approach 
Eoi im. We have then analysed theoretically the effect of time delayed feed¬ 
back from an external mirror on the stability of transverse localized structures in 
a broad area VCSEL. We have derived a real order parameter equation of the Swift- 
Hohenberg type with delay feedback. This analysis is only valid close the nascent 
optical bistability and close to large wavelength pattern forming regime. In this dou¬ 
ble limit, we have estimated the threshold associated with the drift-instability lead¬ 
ing to the spontaneous motion of LS. Explicitly the threshold is given by a simple 
formula Tj T = 1. We conclude that the condition under which transition from mo¬ 
tionless LS to a moving one, depends only on the delay feedback t] and the delay 
time T and not on the other dynamical parameters of the VCSEL system. This con¬ 
clusion is valid for the variational and nonvariational Swift-Hohenberg equations, 
and in reaction-diffusion systems. However, in optics, the role of the phase is im¬ 
portant because the intensity and the phase of the light operate on the same time 
scale. We have investigated in the last part of this chapter the role of the phase of 
the feedback on the mobility of LS. We have shown that, depending on the phase 
of the feedback, it can have either stabilizing or destabilizing effect on the LSs. In 
particular, when the interference between the LS field and the feedback field is de¬ 
structive, the LS can be destabilized via a pitchfork bifurcation, where a branch of 
uniformly moving LS bifurcates from the stationary one. We have calculated ana¬ 
lytically the threshold value of the feedback rate corresponding to this bifurcation 
and demonstrated that the faster the carrier relaxation rate in the semiconductor 
medium, the lower the threshold of the spontaneous drift instability induced by the 
feedback. This is a generic and robust destabilization mechanism in one and two 
spatial dimensions settings and could be applied to a large class of optical systems 
under time-delay control. 

We have described spatially localized structures, recent investigations have shown 
that temporal localized structures have been found in fibre resonators 0104ll90l|l051 
|106l|107lll08l and in VCSELS with a saturable absorber B109I . On the other hand, 
it has been shown that the combinated influence of diffraction and chromatic dis¬ 
persion leads to the formation of three dimensional localized structures often called 
light bullet B110II111II112II113II . We plan in the future to investigate the effect of 
delayed feedback on the spontaneous motion of three dimensional light bullets. 

In order to check our theoretical predictions, we also plan to investigate experi¬ 
mentally the formation of moving LS in VCSELs. In addition, we will analyse the 
role of local polarization dynamics in the formation of LSs in the transverse plane 
of the VCSEL. This would allow us to study the spontaneous motion of vector LSs 
with different polarizations under the effect of delayed feedback. 
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